Effective Quantum Theories for Transport in Inhomogeneous Systems with 

Non-trivial Band Structure 



Christian WickletQ and Wolfgang Belzi^] 
Universitat Konstanz, Fachbereich Physik, 78457 Konstanz, Germany 
(Dated: September 25, 2012) 

Starting from a general iV-band Hamiltonian with weak spatial and temporal variations, we derive 
a low energy effective theory for transport within one or several overlapping bands. To this end, we 
use the Wigner representation that allows us to systematically construct the unitary transformation 
that brings the Hamiltonian into band-diagonal form. We address the issue of gauge invariance and 
discuss the necessity of using kinetic variables in order to obtain a low energy effective description 
that is consistent with the original theory. Essentially, our analysis is a semiclassical one and 
quantum corrections appear as Berry curvatures in addition to quantities that are related to the 
appearance of persistent currents. We develop a transport framework which is manifestly gauge 
invariant and it is based on a quantum Boltzman formulation along with suitable definitions of 
current density operators such that Liouville's theorem is satisfied. Finally, we incorporate the 
effects of an external electromagnetic field into our theory. 
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I. INTRODUCTION 

When performing a semiclassical analysis, one natu- 
rally encounters Berry phased and meanwhile, the im- 
portance of these so-called geometrical phases in con- 
densed matter physics is beyond question^. For exam- 
ple, not long ago it was realized that the electric polariz- 
ability can be defined in terms of a Berry curvature, for 
the first time a consistent formulation of this subjectP. 
Furthermore, research in the field of the anomalous Hall 
effect (AHE) has shown that the intrinsic contribution 
is related to a Berry curvature, which is a quantum me- 
chanical property of a perfect crystaPP Also, magnetic 
monopoles appearing in the definition of a momentum 
space effective magnetic field give important modifica- 
tions to universal conductance fluctuations^. Finally, 
topological interference effects arise from spin Berry 
phases in single molecular magnets^. 

There arc different ways to obtain semiclassical trans- 
port equations (see Refs . EJ an d Bland references therein), 
like wave-packet analysis^^, or the systematic diagonal- 
ization method developed by Gosselin and coworkers^HlSl 
which, however, does not include the possibility of time- 
dependent perturbations. Furthermore, there are vari- 
ous works treating semiclassical quantum transport equa- 
tions that incorporate Berry phase effects: the scenario 
of a general 2-band model is considered by Wong and 
Tserkovnyal^, and spin-orbit coupled systems^ as well 
as a non-Abelean gauge-field formulation^ is investi- 
gated. The Lowdin partitioning, or quasi-degenerate per- 
turbation theory, used in the book of Winkler^ to de- 
rive effective models for certain bands in spin-orbit cou- 
pled semiconductors can also be related to a semiclassical 
treatment, however, there focus is only put on the Hamil- 
tonian, not on the dynamical variables or other aspects 
of the system. 

In this work, we present a self-contained derivation of 
the semiclassical dynamics which is based on the Wigner 



representation^- r phase-space representation of quan- 
tum mechanics - which is a natural starting point for a 
semiclassical analysis. One big advantage is that one can 
obtain corrections systematically to arbitrary order in 
h. Also, a re-quantization of the effective theory is not 
necessary, which is a big drawback of the wave-packet 
analysis which derives a Lagrangian from the equations 
of motion for the wave-packet center of mass coordinates, 
and it is not always clear what the canonical conjugate 
variables are. The relation between canonical and kinetic 
pairs of conjugate variables, however, emerges naturally 
in our formalism. We adopt a 4-component vector nota- 
tion which allows us to incorporate spatial inhomogene- 
ity as well as temporal variation on an equal basis. In 
the course of our treatment, we will find how fictitious 
electric and magnetic fields (real space and its momen- 
tum space pendants) appear in effective theories and we 
complete our work by developing a low energy effective 
quantum transport theory which is manifestly gauge in- 
variant and consistent with a description in the original 
frame. Finally, we address the interesting question of how 
an external electromagnetic field modifies the formalism. 

We have several scenarios in mind of applying our for- 
malism to: studying the electron dynamics in the pres- 
ence of an arbitrary inhomogeneous and time-dependent 
ferromagnetic exchange field, which exhibits many inter- 
esting phenomend2iH23] One can make various general- 
izations like adding a spin-orbit coupling term, which 
would give rise to an even much broader range of new 
effects such as the anomalous Hall effect, and the lat- 
ter would be additionally modified by the inhomoge- 
neous magnetization. Studying spin-transport phenom- 
ena and/or adding thermal gradients delivers a whole 
new range of possibilities.^ Furthermore, there is a re- 
cently discovered class of materials called topological in- 
sulators, in which the quasi-particle momentum is inti- 
mately coupled to the spin-degree of freedom and thus 
behaves similar to a relativistic Dirac Fermion! 25 * 26 ! Now, 
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adding an additional coupling between spatial and spin 
degrees-of-freedom leads to a p hy sically rich system al- 
ready subject to various studies .EzHH] Consid ering the ef- 
fects of mechanical rotation of these system4221Ml j s an _ 
other interesting application for our formalism. 

Approaches that diagonalize band-space in order to 
obtain an effective low energy description have been 
performed in the case of quasi-free electrons in ferro- 
magnetic metals with spatiotemporally varying exchange 
field. There, one can come up with a position-dependent 
unitary matrix U(r) that diagonalizes spin-space which 
maps the problem on that of a free particle in an (fic- 
titious) electromagnetic field.^ The converse situation 
of coupling between spin and momentum as given by 
the spin-orbit interaction can equally be treated by a 
momentum-dependent unitary matrix U(p), for example 
the Foldy-Wouthuysen transformation in the case of the 
relativistic Dirac equation^. 

It is clear that the full quantum mechanical problem 
of surface Dirac Fermions coupled to a general inhomo- 
geneous and time-dependent magnetization texture will 
in general be very complex due to the locking of spin and 
momentum as well as coupling between spin and spa- 
tial degrees of freedom. If one wants to diagonalize spin 
space, the unitary transformation has to involve the pair 
of canonical operators r and p which is rather difficult 
due to non-commutation of r and p. Our approach relies 
on the fact that, if r and p are classical variables, such 
a unitary transformation is much more simple to find, 
and then, we resort to quantum corrections which are 
of the order h. Despite the fact that we are performing 
formally an expansion in h 7 it does not necessarily need 
to be restricted to the semiclassical regime. In fact, as 
pointed out later, our actual expansion parameter might 
be a different one, depending on the physical system and 
the regime under investigation. For example, in the case 
of the Dirac theory, as we will discuss thoroughly in sec- 
tion |VI| the actual scale relevant for our expansion is 
the Compton wavelength A c = — so the resulting Pauli- 
Equation still correctly describes the quantum regime for 
scalar potentials smooth on the scale of A c . 

The outline of this work is as follows: In section [TTJ we 
will introduce the unitary transformation that performs 
a rotation within band-space such that the Hamiltonian 
becomes band-diagonal. Since this transformation is not 
uniquely defined, we will discuss the implications of this 
additional gauge degree of freedom. This motivates a 
description in terms of kinetic variables, which leads to 
the appearance of Berry curvatures, which is discussed 
in some detail in section [nil We will also investigate how 
observables change in the course of the diagonalization 
and we discuss the electronic spectrum as well as energy 
corrections appearing therein. Subsequently, in section 
|IV| we develop a manifest gauge invariant description of 
the physical system restricted to a certain band, i.e. we 
seek a projected theory without the necessity to refer 
back to the original Hamiltonian. To this end, we find 
equations of motion for the quasi-probability density, es- 
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FIG. 1. Illustration of the band-diagonalization scheme: 
Bloch-sphere representation for a simple 2-band model H = 
£(r, p, t) ■ 9 — £ x a x + £ y a y + £ z a z where the vector is rotated 
by Uo(r, p, t) such that it aligns along the z-axis at each point 
in phase-space. 



sentially a quantum Boltzmann equation applicable to 
a non-equilibrium scenario. We also find current densi- 
ties that obey a conservation law corresponding to Li- 
ouville's theorem in classical mechanics. Our quantum 
mechanical equations of motions are formally similar to 
the equations of motion for the center-of-mass motion 
of a wave-packet and, in fact, the latter is just a spe- 
cial case of our formulation. Finally, in section [VJ we 
illustrate how to treat the external electromagnetic field 
which can be done in the spirit of a hierarchy of effective 
theories. After concluding the discussion of general sys- 
tems, in section [VT] we apply the apparatus developed to 
the Dirac equation and readily find a relativistic version 
of the Pauli-Hamiltonian, thereby gaining some interest- 
ing insights into the structure of the Dirac equation. For 
a concise summary of this work, the reader might want 



to go to section VI where all central results are being 
referenced and find immediate application. 



II. THE QUEST FOR A BAND-DIAGONAL 
HAMILTONIAN 



We now consider a Hamiltonian that consists of N 
bands and which is almost diagonal in momentum space 
p, but has some additional spatial- and/or temporal vari- 
ation imprinted on it. In the usual quantum representa- 
tion, it is expressed in terms of the canonical pair of op- 
erators [ri,Pj] = iHSij and additionally carries the N x N 
dimensional matrix structure. Prominent examples that 
fall into this category are the Dirac equation to be studied 
in more detail in section lVTl as well as the aforementioned 
system of surface Dirac fermions coupled to a spatially 
dependent magnetization. 

For practical reasons, we resort to a description in 
terms of single particle Greens functions, or more specif- 
ically the inverse thereof, 5 = ihdtlN — H(r,p,t). One 
reason is that only S represents the complete equation of 
motion, i.e. all kinetic equations involve this operator, 
and not H. alone, for example, we can generally write 
quantum kinetic equations in the compact form 



[S, V] = , 



(1) 
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where T> represents any dynamical variable like the usual 
retarded or Keldysh Greens functions, density of states 
or the density matrix for which ([!]) reduces to the well- 
known von Neumann equation 

iW tP =[H,p] . (2) 

The properties of V (the type of Greens function, etc.) 
enter through appropriate boundary conditions in our 
parameter space. For the example just given, we need to 
fix our density matrix at some initial time p(to) = pq. 

We now transform everything into the Wigner 
representatiorPfil, so that our Hamiltonian and our ob- 
servables are a function of the variables r,p, t and energy 
e so that (see Appendix [A] for more details) 

E(r,t,p,e) = el N - U(r,t,p) . (3) 

Multiplication of operators is now performed by virtue of 
the Moyal product!^, 

it 1 A fA\ 

* = ex P y A i ( 4 ) 

where the differential operator A is given by 

A = %~3 P - d t d e - dpdr + d e d t = - dad* + d^d^. 

(5) 

and we introduced the compact 4-component vector no- 
tation x = (t,r) and 7r = (e, — p). Note that we always 
use x, d x = (dt,d r ) and A x — {At,A r ) in contravariant 
notation whereas we implicitly assume covariant notation 
for the symbols 7r, d„ — (<9 e , — d p ) and A n = (A t , — A p ). 
Since contraction will always be between pairs of x and 
7r, there is no need to indicate covariant vectors. This 
provides us a symmetric and compact notation in the fol- 
lowing treatment. The kinetic equation ([I]) is straight- 
forwardly transformed into the Wigner picture: 

[E* V}=0. (6) 

We are now looking for a unitary matrix U(x, n) that 
transforms our initial Hamiltonian E(x,ir) into a band- 
diagonal Hamiltonian, denoted by S(a;,7r) in the follow- 
ing, i.e. 

U*rf=UUU = I (7) 
U*E*lfl=E. (8) 

Let us note that upon this transformation, U * e * W 
can acquire off-diagonal elements when IA depends on 
time, thus requiring the diagonalization of the combi- 
nation E = el m — T-L rather than H alone. This is one 
major difference to previous semiclassical schemes^. For 
diagonalizations that do not require explicit dependency 
on the time parameter, treating % and E is equivalent, 
and in the following, we will use the term Hamiltonian 
equally for both objects. Note that another advantage of 
using 3 instead of "H, which is the form invariance of the 
kinetic equation (T6J) under unitary transformations. 
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FIG. 2. Illustration of the bands described by the Hamilto- 
nian H and A is a typical energy scale for interband distances. 
The projection operator on the i-th band is denoted by Vi- 
For low-energy processes around the Fermi energy Ef, only 
one or few overlapping bands are relevant and thus contribute 
to, e.g transport properties. Corrections due to the influence 
of other bands enter as Berry curvatures and are oc 1/A. 

In the classical limit, the operators fi and pj commute, 
while in the Wigner representation, the Moyal product 
Q becomes trivial as h — > 0; after all, in this formula- 
tion the Moyal product encodes the non-commutativity 
of the canonical variables. Then, we essentially have to 
diagonalize a N x N matrix "H, whose elements are func- 
tions of r, t and p. We call the unitary matrix associated 
with this rotation in band-space Uo(r,p,t), so that 

U EUl = S «. u nul = n , (9) 

and the digonal elements of So constitute the classical 
energies of the N bands described by our initial Hamil- 
tonian. At any rate, we assume from now on that we 
know the diagonalization matrix Uq analytically. Note 
that all our matrices parametrically depend on x and 
7r, so we essentially diagonalize locally at every point in 
2 x (3 + l)-dimensional parameter space (see Fig. [ljfor an 
illustration), which becomes meaningful in the semiclas- 
sical limit, as position r, momentum p, time t and energy 
e are well-defined in this limit. Also note that Uq (r,p, t) 
is not uniquely defined which brings up the problem of 
gauge-invariance as discussed in detail later. 

Furthermore, if the band degrees-of-freedom couple to 
only either x or 7r, then we already found the exact ex- 
pression U = Uq. For example, treating metallic ferro- 
magnets with inhomogeneous magnetization and neglect- 
ing spin-orbit interactions^!, it just necessary to only di- 
agonalize m(r)<r, since the energy dispersion is diagonal 
in spin-space. Then, the effect of IA is that p acquires an 
additional vector potential: p — > p — iftlA(d r W). 

The general situation, where the degrees-of-freedom 
couple to both spatial and momentum coordinates si- 
multaneously, is much more involved and chances that 
we can come up with an exact solution for U(x,tt) are 
slim. Therefore, we adopt a gradient expansion approach 
where we expand the Moyal Product * in powers of h 
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and calculate corrections systematically order by order. 
The fact that we are performing formally an expansion 
in h does not necessarily imply that we are restricted to 
the semiclassical regime. Our actual expansion parame- 
ter might be a different one, depending on the physical 
system and the regime under investigation. Typically, 
the expansion parameter is the ratio of the energy as- 
sociated with band-dependent spatiotemporal variations 
and the interband energy distance A, as illustrated in 
Fig. [2] and therefore, the larger the band separations (or 
weaker coupling), the better the approximation becomes. 
To emphasize this fact, we will introduce the notation of 
pc in this article to indicate that the expansion is not 
necessarily a semiclassical one. 

We assume that we already solved the zeroth order 
problem which yields the unitary matrix Uq, so that now, 
we can introduce gradient corrections to this matrix by 
writing 



Hcrmitian 



U=(l +U X +U 2 



(10) 



which describes the matrix that diagonalizes our Hamil- 
tonian H to arbitrary order in pc. IA\ is oc pc and will be 
specified explicitly later. As a result of the expansion of 
the Moyal product, we obtain the covariant derivatives 



di — UodiUl = di — iA 
which acquire the Berry connections 



(11) 



(12) 



where di is a placeholder for any of the possible deriva- 
tives d t ,d e ,d r and d p . Berry phases, or often called ge- 
ometric phases, are omnipresent in modern physics, and 
describe phases picked up along a trajectory in curved 
geometries. Note that Ai is a N x N matrix which is 



A] 



i{U diUiy = -i{diU Q )ul = iU {diUl) = Ai 



(13) 



where the last step is due to = 9,;(WoWq) = (diUo)Ul + 
Uo(dilAl) and the diagonal elements of Ai describe the 
usual Berry connections arising from transport within a 
certain band. The off-diagonal elements mix contribu- 
tions from two bands and thus describe effects due to 
inter-band coupling during transport in a certain band. 
As we will see later, these inter-band transitions will 
give rise to important corrections for example the cor- 
rection to the energy that appears to leading order in 
pc. Transforming higher order derivatives is straightfor- 
ward, except for the additional freedom of exchanging 
partial derivatives, which leads to a general relation be- 
tween Berry connection matrices A4. Starting from 



UodjdiUl 



and likewise, 



djdi 



(dj - iAj)(di - iAi) 



(14) 



UodidjUl = didj = (di - iAi)(dj - iAj) , (15) 

and, if we assume symmetry of second derivatives djdi = 
didj, the following relation should hold 



diAj 



djAi 



i \Ai , Aj 



(16) 



In fact, this identity can be directly shown by using the 
definition of A ( 12 ) and the exchange of partial deriva- 
tives. 

Now we have everything at hand to systematically cal- 
culate S = So + Si + S2 . . • to arbitrary order in pi and 
to first order, we explicitly obtain the following result for 
the transformed expression, 



Hi — U\Cq + ^qUI — — {Atv, d x ^o} + 77 {A x , <9tt"o} + {A^qAx — Ax^oA 



ih 



(17) 



We obtain an additional constraint for U n from the con- 
dition of unitarity Q , which to first order in pc reads 



Ul +U\ + y [Ak,Ax] = , 



(18) 



with y\ = —y \ assumed to be antihermitian, such that 



condition ( 18|) is satisfied. Plugging U\ back into Eq. 
(fl7| yields 



and is obtained by substituting 1 for ^0 into the trans- 
formation relation Eq. (17). Since this relation fixes the 



Hcrmitian part oiUi, we can make the ansatz 
ih 



[A 7T ,Ax] + yi , 
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li = [y u H ] - \ {A„, [d x - iA x , So] } + ^ - §Ar, S ] } + 0(^ 2 ) 



(19) 



We observe that [3\,Ho] is completely off-diagonal and 
3^1 is well suited to absorb the off-diagonal part of the 
last two terms on the right- hand-side of Eq. (19), in the 
following denoted as 1Z . After all, the objective is to 
diagonalize 2, so we want the off-diagonal part of Si to 
vanish, which is achieved by the condition 



[3?i,Ho] =K . 
The solution can be readily given as 



(20) 



(21) 



in the case of a completely diagonalized (Hq)^ — ti&ij 
and which is antihcrmitian as desired. Here we explicitly 
see that corrections due to y are inversely proportional to 
the separation between bands, and thus are of the same 
order as corrections due to Berry phases. In fact, the 
commutator on the left side with the Hamiltonian Hq is 
exactly what also appears in the von Neumann equation 
suggesting further that this term describes corrections 
due to inter-band dynamics. 

We remark that the difference between the diagonaliza- 
tion of S and H is merely the modification of 3^1 by that it 



acquires an additional contribution from the off-diagonal 
elements of the Berry connection At ■ However, this term 
is crucial when our transformation 1A is time-dependent 
and only then will the formalism yield consistent results. 
Note that the definition ( 20 1 does not fix the imaginary 
diagonal part of but without loss of generality, we 
can set this part to zero. A non-zero imaginary diago- 
nal part corresponds to the linear ord er expansion of the 
gauge phase factors e lx< - x ^ (see (23 1). Essentially, the 



freedom of choice here can be reduced to the problem of 
gauge invariance, to be discussed in the next section. 

Yet, we also want to capture the situation of degen- 
erate or overlapping bands, so the desired form of S is 
in general block-diagonal. To formally express this mat- 
ter, we introduce the projectors Vi that define the bands 
(see Fig. 2| and that, in the end, we want use for our 
effective theory. Of course, the set of Vi has to be spec- 
ified together with Uq, since the unitary transformation 
has a freedom of how we distribute the bands amongst 
the entries of our matrix. For example, a natural choice 
would be to sort the bands with respect to their ener- 
gies. Then, introducing the projected Berry connection 

A^ = VdA^Vd = ^iViAjrVi, we can write for the 
Hamiltonian that is diagonalized up to first order in pt, 



S = S - JU^flUSo 



HA^d n E + jV d ({An, [A X ,E ] } - {A*, [X.Ho]}) V d 



(22) 



As compared to expression (191, we can drop two anti- 
commutators, since after truncation, A^ and So com- 
mute. As discussed more thoroughly below, this is even 
the case when either So is block-diagonal and/or the pro- 
jected Berry connections are non-Abelian. 



Apart from the obvious term So = e — which can 
be understood as the classical energy, the last term of 
(22 1 can be thought of as correction to the energy due 



to inter-band transitions and corresponds to the energy 
associated with persistent circulating currents, as for ex- 
ample the magnetic Zeeman energy in the case of the 
Dirac equation (see discussion in section VI I . The two 



middle terms in ( 22 1 appear to be Berry phase correc- 



tions to this energy, however, they are not unique in the 
sense that they depend on the specific form of Uq. The 
meaning will become more apparent later, but before, we 
address the question of gauge invariance of the effective 
Hamiltonian (22 1. 



A. Gauge Invariance 



We mentioned previously that there is an additional 
degree of freedom in the choice of unitary transformations 
Uq and Uq which all yield the same diagonal Hamiltonian, 
S = WqSWq = UqEUq . These different unitary transfor- 
mations are related by local phase factors in 2 x (3 + 1)- 
dimensional parameter space, so that we can formally 
connect two unitary transformations hi and W by the 
(block-)diagonal N x N phase matrix 



$(a;, it) = 



/ e ixi(a>.T) o 

e *X2(x,7T) Q 

o e *X3(v,Tr) 



V 



(23) 
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where ^„(s, tt) are arbitrary functions in our parameter 
space, which are appropriately termed gauge fields. In 
fact, W = &U describes a gauge transformation and 
can be thought of as a local phase transformation in an 
extended phase space which includes time and energy. 
This gauge transformation changes the Berry connection 
matrices according to 

A' k = iU' Q d k U$ = $A k & + i$d k & = ®A k & + X k , 

(24) 

where the field X k = i$d k & = <9 fe diag(xi, X2, Xn) is 
a (block-)diagonal matrix containing the partial deriva- 
tives of the phases. &Ak& modifies the off-diagonal el- 
ements by giving them additional phase-factors, while 
the change in diagonal elements is due to X k . In prin- 
ciple, we have to distinguish three different cases here, 
the first one being that So is completely diagonalized 
with well-separated bands which corresponds to the sit- 
uation just described. However, when we keep part of So 
block-diagonal because bands are overlapping or degen- 
erate and we can distinguish the bands, i.e. they have a 
physical meaning that we want to retain (for example we 
have spin-degenerate bands but want to describe spin- 



dependent physics) then we have only a U(l) gauge free- 
dom within this sub-block. Third and last, if we have M 
degenerate bands, i.e. there is a M-dimensional sub block 
in So that is proportional to the unit matrix, and further- 
more, we cannot or do not want to distinguish between 
the degenerate bands, we have the additional degree of 
freedom to rotate within this degenerate space giving us 
an additional SU(M) gauge invariance. Contrary to the 
first two cases, this last one describes a situation with 
the effective description of this M-dimensional subspace 
being a non-Abelian gauge theory with the symmetry 
group U(l) x SU(M) and consequently, X k constitutes a 
non-Abelian field. We do not differentiate between these 
cases explicitly in the following because they are straight- 
forwardly treated in our formulas, thus requiring us only 
to comment in situations where special care is required. 

According to the preceding discussion, the projected 
Berry connection matrix is only modified due to X , 

A'^ = V d A' x V d = A^ + X x , (25) 

and the alternative transformation due to Uq leads to the 
Hamiltonian 



= S - \ {A^ d \d x E } + \ {AW&So} + jV d ({At, [A x ,E ]}- {A X1 [A«, S ] }) V d + O(^) . (26) 



The last term, representing the inter-band transition cor- 
rections to the energy, does not change, since the addi- 
tional terms due to X Xl X n are projected out by V d and 
thus, as a consequence, are absorbed by y[ viz. U[ which 
of course does not need to coincide with U\ . 

However, the other two terms linear in do explic- 
itly depend on the gauge, and therefore change the ef- 
fective Hamiltonian. Clearly, this shows that the effec- 
tive Hamiltonian alone is an incomplete description as 
it directly depends on this additional degree of freedom. 
Therefore, in order to make any sense out of this, we 
have to identify our physical observables, because in the 
end, the physical results derived from our effective theory 
should not depend on a specific gauge. 



III. CANONICAL VERSUS KINETIC 
VARIABLES AND GAUGE INVARIANT 
DESCRIPTION 



The previous section showed us that there is still an 
ingredient missing in our effective theories. In order to 
investigate this matter, let us study the dynamics of our 
system, and construct the effective theory such that the 



results obtained within this description are consistent 
with what one would obtain in the original frame. The 
question is now, whether one can find a manifest gauge 
invariant formulation and, as we will explain in the fol- 
lowing, it is indeed possible. 



Parameter transformation to kinetic variables 



Let us consider the operator S that describes some 
physical observable of our system and, in performing 
the rotation that brings our Hamiltonian T-L into diag- 
onal form, it transforms our observable S along with it. 
The observable in the rotated frame S is then related to 
the original operator S by virtue of relation (48 1. We 



now consider the projected system, i.e. we have in mind 
to develop an effective, yet exhaustive description of the 
physics taking place within a certain band that is suffi- 
ciently well separated from all other bands in order to 
treat this band independently to a good approximation. 
For the moment, let us assume that our observable S is 
band-diagonal, i.e. it is a scalar function of x and 7r, so 
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that the last term in relation ( 48 ) vanishes and we have 



s = s - ~ {Ai d \d x s Q } + h - [a£\ d^So] + oo?) 



(27) 



where Sq = UqSUq. Later, we will lift this restriction and 
consider a matrix S with general structure in band space, 
so that we will also get this additional term, giving rise to 
important contributions. However for the moment, (27) 



is nothing but a Taylor expansion of So to first order in 



S = S Q (x- KAg\ir + hA (d) ) + 0(p?) 



(28) 



suggesting the parameter transformation to band pro- 
jected kinetic variables X and II 



= x- hA^ 



so that we can write 



s = s (x,n) + o(p 2 ) . 



(29) 
(30) 



(31) 



We note that to leading order in p, we can equally write 
A(x, 7r) = A(X, H)+0(p), since Berry connection terms 
are already linear in p. 

P = p — HA r is of course well known in the Hamilton 
formulation of particles in an electromagnetic field. In 
an analogous way, the position R = r + hA p acquires an 
additional Berry connection with its corresponding Berry 
curvature, or momentum space magnetic field that gives 
rise to the so called anomalous velocity term. Further- 
more, E = e + HAt attains a contribution which gives rise 
to an electromotive force and appears in the form of an ef- 
fective electric field (see for example Ref. f3"5j) . The same 
is true for the electromagnetic field, where the electric 
field can be also rewritten in terms of a time-dependent 
phase, effectively changing the gauge. Finally, for rea- 
sons of symmetry, one would also have T = t — hA e but, 
at least in non-interacting Hamiltonian systems, A t is 
zero, since the energy-dependency in S is trivial. How- 
ever, the situation is different if one considers an inter- 
acting system and uses an effective non-interacting quasi- 
particle description, because the self-energy that includes 
these interaction in general carries a non-trivial energy- 
dependence inherited by the single-particle greens func- 
tion and thus by Sp^l 

Let us have a look at the commutator relations between 
the kinetic variables) 13 * 15 ! 



R 



(d) 



R 



(d) 



R 



p(d) * p 



,(d) 

j 

(d) 
'3 

(d) 



= iPe ijk Bf 



(32) 



where we introduced 

s M = h ( d r x A {d) 



ih 



(a^ x 4, d > 



(33) 



B (p) = h (d p x A { p d) j - ih [A^ x j , (34) 

which can be considered as a generalized magnetic field 
in real space and reciprocal (or momentum) space. Such 
non-Abelian Berry curvatures have been treated in Ref. 
[TS1 Furthermore, 



(9A, 



(d) 



dA 



(d) 



dp. 



Aid) A(d) 

T"i 1 J) ■ 



(35) 



As we will also encounter later, the dimensionless tensor 
describes the change in the metric of the phase-space 
due to the parameter transformation from canonical to 
kinetic variables ( 29 1 . We can make this more apparent 



by relating it to the change in differentials 



di? ? ; 



hd r j A Pi ) dr 



so that 



dPi = (Sij - hd P] A Ti ) dpj , 



dR dP = dr (1 + 0") dp + 0{p?) 



(36) 
(37) 

(38) 



A more compact way to write these Berry curvatures 
is to use the covariant derivative ( 11 ), projected onto our 
band-diagonal space 3^ = VddxVd, for example 
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(8 { t d) A^ _ BlfAif 



(39) 



It is well known that the Berry curvatures are invariant 
with respect to gauge transformations, and the commu- 
tator is essential as it provides the full SU(M) gauge 
invariance in the non-Abelian case. 



In accordance with the effective magnetic fields (33) 
and (34), we introduced the effective electric fields 



(d r A[ d) -d t A^-i[A^\4 d) }) 



£ (p) = h (d p Af ] - d t A {d) - % \A {d \A 



(40) 
(41) 



which shows us indeed that At appears in the role of a 
generalized electric potential, however, it can also depend 
on momentum P. In the Abelian case (for non-Abelian 
fields, it works if we take Trj\,/£ and Tr^S or when we 



take the covariant derivatives ( 11 ) along with the full ma- 



trix structure of the Berry connections) , these fictitious 
fields obey homogeneous Maxwell equations 



a 



r/p 



d r/p -B (rM =0 



(42) 



however, in order to determine these fields independently 
as in classical electrodynamics, one would need two ad- 
ditional inhomogeneous equations containing (effective) 
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source terms as inhomogeneities. Note that in general, 
2? (7,/p> and £ (rlv) depend on r and p simultaneously. As 



in (32 1, the effective magnetic fields can be also defined 



in terms of commutator relations 

£M _ J_ Jtj;(d) * p(d) 

ih L 

gWi = _1_ \ E (d) » H (d) 



(43) 



There exists a sum rule for the fictitious fields 



Tr N B lp) = 
Trjv£ (p) = 



c N B M = , TrAre^ = , 
?n£ m = , (44) 



which can be found by taking the trace over all bands (we 
denote this sum over bands as Trjv here and throughout 
this work) and making use of the identity (16). This 
means that all the effective forces for each band balance 
each other in total, or in other words, if all bands are 
completely filled (and thus, the density operator p is pro- 
portional to the unit matrix In), the system does not 
experience any net force. 

The basis of all calculations within the Wigner repre- 
sentation of quantum theory is the Moyal bracket \_S * T] 



between two operators S and T . Now, we want to rewrite 
this in terms of kinetic variables only, which is achieved 
by transforming the derivatives in the Moyal product to 
act on kinetic variables. For concise notation, we in- 
troduce the tensor of Berry curvatures as (a, j3 denote 
indices with respect to (t,r,p)) 



'A w A (d) 



(45) 



or explicitly in terms of the fictitious fields used previ- 
ously 



e 



-£<"> 

-(er-) 1 




(46) 



e - i R <p > 



where we recognize the top-left part of O as being essen- 
tially the electromagnetic field tensor. 

Neglecting terms of order 0(pc 2 ), we find the explicit 
form of the Moyal product after the transformation to 
kinetic variables 



exp 




(47) 



This is a central result of this section, since it shows that 
expressing all quantities in terms of kinetic variables al- 
lows us to deal solely with manifest gauge invariant ex- 
pressions. Essentially, the bottom line of this parameter 
transformation is that it changes the metric of the Moyal 
product by the appearance of the Berry curvatures 0. 
Immediate consequences of result ( 47 ) are the equations 



of motions to be discussed in detail in section [TV 



B. Transformation of general operators 

Let us briefly comment on observables with non-trivial 
band matrix structure 0(x, 7r), and their transformation 
into the rotated frame, which is performed analogously 
to the transformation (19), 



o = u * o * = 6 + [y x , 6 ] - 1 {A,,, [flu - f A*, do] } + \ {a x , [a* - o ] } + o(p?) , (48) 

and the back transformation is given by 

Uooul =u (uU6* u)i4 = 6-h[y 1 ,d] + ~{A 7r , [a m -%A>,6\}-\ {a x , [d„ -^,o]} + o{n 2 ) , (49) 



r 



which can be readily checked by plugging Eqn (48) into pL . Now, treating the diagonal and off-diagonal part of 



the back transformation ( 49 ) and dropping terms of order O = UqOUq = O q + O separately, we find 
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1 } ~ V 



4 (5 (d) 



}]^ d + 0(/r 2 ) , 



(50) 



and the contribution arising from the off-diagonal part of Oq , 



h 

2 



{ A \b _ i/l I + 5 / A 

y *^7v i u x 2 35 5 J O |^ ' 



<^7T 2 •^ 7r ' 



}) V A + O(^) 



(51) 



which are both independently gauge invariant. While the 
gauge invariance of the former is straightforward to show, 
the later requires more work and we have to take into 
account that y± is modified under a gauge transformation 
as 

yi -> yi - ~ ({xxA} - be*, a*}) , (52) 

along with 9 Q (5 — > <9q,C> + i [x a ,0] and _4 Q — > ^4^. + x Q . 



C. Expectation values in the rotated frame 

Let us now go back to the initial question of the dynam- 
ics of our system within the effective theory by studying 
the expectation values of physical observables, which can 
be obtained in the Wigner representation by the integra- 
tion over the complete phase space, 



(S) =Jd d rJ Trjv {p(x, tv) * S(x, ir)} , 



(53) 



and the trace is with respect to the matrix structure. 
Note that the factor l/(27i7i) d = l/h d describes the 
proper quantization of the phase space volume and thus 
is directly obtained by transforming quantum averages 
into the Wigner representation. 

If we assume the integration over the whole phase space 
to be unbounded and any surface contribution from the 
integrand at infinity to vanish, we can perform partial 
integrations to show that all the partial derivatives in 
the Moyal product * cancel each other, so that we can 
equally write 



(S) 



d d r 



d d p 



Tr N {p(x,Tr)S(x,ir)} . (54) 



By using the cyclic property of the trace and by partial 
integration we can easily show that 



(O) 



{2-kH) 

and we can write 
(S) = (S) = [d 



' [ldd ^Tv N {U*0*U^} = (6), (55) 



d d p 
(2ttH) 



- d Tr N { P (x,tv)S(x,tv)} 



(56) 



In equilibrium, the density operator is given by 

p(x, tv) = / D (e) *(H) - / D (e) 5(e - H(x, tt)) , (57) 

or, if we do not want our results to be energy resolved, 
we use directly 



p(x,p)= J de f D (e)5(e-H(x,7r)) 



(58) 



which gives us the density in phase space and, as we 
will discuss later however, it is to be interpreted as a 
quasi-probability. An explicitly time-de pen dent Hamil- 
tonian H(x,Tr) has to be treated using (96) instead. If 
a band is completely filled, it becomes 1 at the diagonal 
element corresponding to that band. If energies are well 
separated, and we assume excitations localized in energy 
space, we can assume any off-diagonal elements in the 
density operator to vanish. In fact, those off-diagonal 
entries correspond to coherent excitations that are split 
amongst several bands. If the band splitting is suffi- 
ciently large, this leads to rapid oscillations subject to 
decoherence. In the end, the projection operation defined 
by Vd is essentially enforced by the diagonal representa- 
tion of the Hamiltonian S, and by the density matrix 
p which gives us only those states that have coherences 
within bands (or degenerate/overlapping bands so that, 
again, the energy argument applies). In particular, this 
is certainly true for low energy transport, where physics 
takes place only in the vicinity of the Fermi level. 

In order to proceed, let us transform the integration 
variables to kinetic ones, and in doing so, we also have 
to take into account how the volume element in phase 
space changes, which is given by the determinant of the 
Jacobian 



D~ 



dot 



f 



d(R,P) 



TrG 7 " + 



det (1 + e^) 



(59) 



or D(X, II) = 1 - Trg^ + 0(pt 2 ). In non-Abelian situa- 
tions the non-trivial matrix structure of the Berry curva- 
ture will be inherited by D which will be accounted for by 
performing the integration before taking the trace, thus 
yielding 



(S) = Tr 



JY 



d d R d d P 
(2nH) d 



D(X,U) p(X,U)S(X,U) 



(60) 
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D describes for example charge accumulation in the case 
of a topological insulator with a magnetization structure 
induced by ferromagnetic exchange. This effect of the 
Berry curvature D on the density of states has been al- 
ready discovered by Xiao and coworkers^. 

Let us look at this in another way by using result 



density (see Eq. 84). With a series of straightforward 



( 49 1 to transform an observable O back into the original 



frame. In addition, we let S be a general observable that 
can posses an arbitrary matrix structure, so that contrary 
to relation (27), the additional dipole term becomes rel- 
evant. In the end, we want to establish the connection 



with (60), so we are interested in expectation values or 



phase-space densities (which then are quasi-probability 
distributions as discussed later), 



s(x, tt) = Trjv I I { P * S} ) = TrjvO 



(61) 



We identify O = \ {p * S} and plug in the back trans- 
formation ( 49 ) , so that 



s(x, tt) = Ty n {O + hAAO - hA^O) 

+ ihTr{6[A x ,A jr ]} + 0{n 2 ) . 

According to the discussion above, it is reasonable to as- 
sume that our observable is given as a function of the 
kinetic variables, i.e. 0(X,T1) and it is instructive to 
treat the inter-band and the intra-band contributions 
separately by splitting 6 = V d OV A + 6^ = 6^ + 
and likewise for S — ^-S^ .so that the contribution 
from the diagonal part becomes^ 



,(d) 



n 



0{p?) , (62) 



which basically undoes the variable transformation so 
that we go back to the canonical pair of variables and 
can write 

s^(x,tt) = Tr N {p(x,ir) D(x,tt) S {d) (x, tt)} + 0{U 2 ) ■ 

(63) 

In addition, we rewrote the last term with the help of 
identity (16), i [Ax^A-k] — d x A^ — d^A x and, according 



to our previous discussion, we have 9 £t = 0, as was as- 
sumed to not explicitly depend on the energy parameter 
so we can replace this term with Tr0 rp . This contribu- 
tion gives rise to the correction factor D(x,tz) that we 
already encountered before, and thus, the last result is 
consistent with relation ( [60| . 

The implications of the diagonal part of S can be 
summarized as undoing the parameter transformation 
together with the appearance of the correction factor 
D(x,p) which locally changes the density. However, it is 
not always possible to ignore the off-diagonal part of the 
observable S, one prominent example will be the current 



manipulations involving the cyclic property of the trace 
along with identity (16), we eventually arrive at 



a^(x,ir) = d a Tr N p-{A„,S^} 



d„Tr N p-{A X ,S^} + 0{t?) , (64) 



so that both contributions to the expectation value, (63) 
and (64) together read 



-d n Tv N p-[A Xl S^}+0(fr 2 ) . (65) 

The importance of these last two terms will become 
clearer in section IV A| when discussing the kinetic equa- 
tions of the effective theory. 

To summarize this section, we have seen that when 
we use the kinetic terms X and II as basic quantities 
for our observables, we end up with expressions that are 
manifest gauge invariant (c.f. Eqns (47) and (60)). In 



fact, these kinetic variables appear consistently in virtu- 
ally all equations of physical relevance, and moreover, it 
is exactly these quantities that we obtain, if we transform 
the canonical variables into the rotated frame, 

X =U * x *W f = x - ihU * dJA^ 

n = u * tt *w f = tv + %m * d x u^ . 



D. Is the diagonalization transformation canonical? 

Before continuing, we would like to point out that with- 
out the projection, X and II still obey the canonical 
commutation relations, which can be easily seen by not- 
ing that the set of unprojected Berry curvatures vanishes 
according to identity ( [l6| (for example, (34) or ( [41] ) with 
^4( d ) replaced by the N x N Berry connection A) . This is 
actually not surprising, since then our unitary matrix Uq 
is connected by SU(N) gauge invariance to the identity 
transformation which clearly has vanishing Berry curva- 
ture. Or in other words, U itself is a SU(7V) gauge- 
transformation which naturally keeps the full SU(7V) 
Berry curvature invariant. This implies that our uni- 
tary transformation is a canonical one, however, since X 
and II now posses a complicated matrix structure in the 
iV-dimensional band space, they no longer commute with 
non-trivial matrices within this band space. For example, 
in a simple particle-hole symmetric two-band model 



T-L = E(p) ■ a + V{r)t 2 , (66) 
our band-diagonalized Hamiltonian has the form 



■H = E{p)<r z + V(R)l 2 , 



(67) 
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while R — r + hA p acquires a 2 x 2 Berry connection The explicit structure of fi is 
matrix, so that [Ri,Rj] = = [Pi,Pj] and [Ri,Pj] 



ihSij. Instead, the commutator [R, a z ] = H[Ap,a z ] ^ 
now encodes the complicated dynamics of inter-band 
scattering, making the problem as a whole not easier 
tractable so, in the general case, the only way out is the 
truncation scheme. And only due to the restriction of 
the Berry connection matrices into a certain sub-space do 
the corresponding Berry curvatures yield a non- vanishing 
value. 



fi 




(70) 



which has been defined in analogy to the Berry curvatures 
given in ( 46 1 . 

Now, we can express the energy term as 



H M = Trfi 7 " , 



(71) 



E. Electronic Spectrum and Magnetic Dipole 
Energy 

In course of the preceding discussion, we have seen 
that rewriting the Hamiltonian in terms of kinetic vari- 
ables would render it gauge invariant. However, in or- 
der to calculate the electronic spectrum, one essentially 
has to fully diagonalize it, and which has to be done in 
terms of canonical variables. But it turns out that this 
remaining gauge-dependence of the Hamiltonian in the 
canonical representation would only affect the wavefunc- 
tions or quantities that build upon them like the retarded 
Greens function or the density operator. These objects 
will acquire local phase-factors that depend on x and p 
and we will illustrate this point later by explicitly study- 
ing the situation in the case of the Dirac equation. At 
this time, we can conclude that the electronic spectrum 
of the system is also gauge-independent. 

Let us briefly discuss the last term of the transforma- 



tion equations like (22), and that we ignored by now, 



n M = jV d ({A r , [A p ,Ho] - {A p , [Ar,Ho\}})V d 



(68) 



which gives rise to corrections due to virtual transitions 
to other bands. The term "magnetic" comes from the 
fact that Hm corresponds to the energy of a magnetic 
dipole in an external magnetic field, where the magnetic 
dipole is an intrinsic property of the band. For exam- 
ple, in the case of the Dirac equation, this term becomes 
essentially the magnetic Zeeman term, though in other 
scenarios, one obtains generalizations thereof, and even 
in the absence of external magnetic fields, JIm can be 
non-zero. 

Since we will encounter terms like in Hm later, let 
us define the quantities (a, (3 are any combinations of 
(t,r,p)) 



a, 



ih 



V d {A a , [Ap,Ho]}V d , 



(69) 



which is an antisymmetric tensor in (t, r,p)-space, and in 
addition, it is a band-diagonal matrix, or block diagonal 
according to the structure defined by the projector Vd- 



where the trace is only with respect to coordinates, and 
not band indices, and we write explicitly 



a- 



-V d {A ri , [Ap^Ua^Vd 



(72) 



As will become more apparent later, terms involving the 
fi tensor are related to circular currents and give rise to 
important terms that should not be ignored. 

For the sake of completeness, and in order to estab- 
lish the link to other treatments in literature (for ex- 
ample Ref. [3] and references therein), let us express 
our quantities and fi in terms of Bloch functions 



icL 



with the Bloch band indices i and 



j. We start by introducing the gauge-invariant transition 
elements 



r 



= -2m{(,Aa)ij{Afi)ji} 

-2h% {(ui\id a \uj) (v,j\idp\ui)} 



(73) 



which are anti-Hcrmitian, T 



= -T { ^ f3) due to the 



Hermiticity of A and, as the name suggests, they describe 
corrections due to virtual transitions between band i and 
j. Then the Berry curvatures projected onto band i can 
be expressed as 



h{d a Af 



(74) 



and one can interpret the Berry connection as being cor- 
rections to the kinetic variables x, tt due to virtual tran- 
sitions into all other bands j . The sum rule ( 44 ) is then 



a direct consequence of the anti-Hermiticity of T. In the 
same spirit, the dipole terms projected on Bloch band i 
read 



fi 



a/3 



IT 



(«,/; 



(75) 



where the energies are = (Ko)«- 

As might be apparent from the above expressions, 
there exists a direct link between and fi in the case of a 
two-level system, and which we establish in the following. 
First, a general diagonal 2-level Hamiltonian Hq can be 
divided into a symmetric part cx I2 which consequently 



commutates out in expression (69), and an antisymmet- 
ric part cx a z responsible for a finite contribution. After 
some simple algebra, we find Tr2(T 2 fi = which implies 
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that is band diagonal, and the two bands get the same 
contribution, in particular the magnetic dipole energy 
is the same for the two bands. The diagonal part is ob- 
tained by taking the trace and using its cyclic properties, 



IV. EQUATIONS OF MOTION 

We are now interested in the dynamics of the sys- 
tems, especially in a consistent effective descriptions of 
the physics within a certain band. 



= hTv 2 H (d a Af ] - d p A$) = Tr 2 (H 9a,p) , 

(76) 



where in the last two steps, we used the identity ( 16 ) 



along with definition ([74]) to replace the commutator. We 
can finally cast the result into the form 



0=2 Tr2 (^°i) 12 



(77) 



where in this expression, the diagonal contribution of H.Q 
drops out because of the sum rule Ti'20 = 0. 

In particular, we can express the magnetic dipole en- 
ergy as 



n M = 2 Tr 2 (£ Tre*") 1 2 



(78) 



where TrB^ = Y?,i®u> whereas Tr 2 denotes the trace 
with respect to the two level band space. 



A. Quasi-probability distributions and currents 

Within the Wigner framework, the basic quantity is 
the quasi-probability distribution p(r,p), while the usual 
momentum and position probability distributions can be 
generally defined as marginals 



p{r) 
P{P) 



d d p 
{2-kK) 
d d r 

(2ttH) 



p(p, r) 
d P(p,r) 



and they are linked to the probability interpretation of 
quantum mechanics, i.e. p{r) — |'0( r )| an d p{p) — 
\ip(p)\ 2 - In the rotated frame however, one has to ex- 
press them in terms of kinetic variables instead, in order 
to obtain gauge invariant results. 

Now the goal which we want pursue in the following is 
to find proper quantities in the rotated frame that lead 
to a consistent physical description when projected onto 
a certain band. 

Since the velocity operator is given by 



dR 

~dt 



1 



[§**] 



and inspired by result ( 65 ) for the expectation values in 



the rotated frame, we define the quasi-probability density 
n, current density j and force density q as follows 



n(R, P,t) = Trjv {p D} , 

j(R,P,t) = Tr N jp - V R xTT N p^ - V P -Tr N prr , 



q(R,P,t) = Tr N ip D 



.dP 
~dt 



V P x TrAr pCV* + V R ■ Tr N p(Q T 



(79) 
(80) 

(81) 



where we used that = — (f2 77 ') T . These densities obey 
a conservation law in phase space, or Liouville's theorem 
that states 



d t n(R, P, t) + V R j{R, P. t) + V P g(R, P,t) = 



(82) 



We substitute the definitions ( 79 )- ( 81 ) into Liouville's 



theorem and establish the identity ( 82 ) after some alge- 



bra, using the kinetic equation ( 94 ) for p and the equality 



of second partial derivatives and dropping terms of order 



0{h 2 ). Furthermore, we need to use the following iden- 
tities between Berry curvatures 



d t D - V R S ( ^ + V P £ M = 



-d Pi D-d Pk QZ-{V R xB 







(83) 



which can be readily shown by plugging in the defini- 
tions Q-d35}, Q and plj). This result strongly em- 



phasizes the importance of including the correction factor 
D(R, P, t) into the expectation values. 
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Here, it is important to realize that the current matri- 
ces have an off-diagonal structure that one has to take 
into account in the light of expression ( |65[ ) and which 
eventually leads to the last two terms in the current den- 
sities (80) and (81), and constitute divergence- free, or 



circular currents. Explicitly for the 4-component current, 
this off-diagonal part takes the form 



dX\ 



[T-Lq , A T , 



(84) 



and is completely expressed in terms of quantities of the 
rotated frame and arises due to the off-diagonal elements 
of the Berry connection matrix. 

Just like the Wigner function is a quasi-probability and 
can be given physical sense only after taking expecta- 
tion values, the same applies to n, j and q. In partic- 
ular, physical meaning can be given only to quantities 
like n(R) or j(R). This is related to Heisenberg's un- 
certainty which states that momentum and position un- 
certainty have to be larger than Plank's constant h, i.e. 
APAR > h. The same is true for the conjugate variables 



time and energy. 

Liouville's theorem constitutes a conservation law for 
the quasi-probability densities, and one can for example 



integrate ( 82 ) over all momenta in order to find a conti- 



nuity equation for the probability densities of the particle 
and current densities, 



d t n(R,t) + V R j(R,t) = -<f d d - x P q(R,P,t) 

JdV 



(85) 



where the closed surface integral on the right-hand side 
has been obtained by virtue of Gauss's theorem, and it 
vanishes if we assume that there is no net momentum 
current flow through the surface at infinity. Explicitly, 
the particle density is 



(2nH) 



n(R,t) -- 

and likewise for the current density. 



d Trjy (p D) 



(86) 



j(R,t) 



ddp t f-n dR 



dt 



Vrx 



d dp 



(2irh)° 



Tr npQ t 



(87) 



while the previous discussion has shown that these defini- 
tions are in fact meaningful physical quantities and con- 
stitute the central result of this section. We will see in- 
teresting implications of the additional dipole term in the 
current when studying the examples of the Dirac equa- 



tion. 



For the sake of completeness, we use result (65) to 



define densities corresponding to the physical quantity S 
in the diagonalized frame 



S(R,t) EE 



d d p 
{2irhY 



Tr 



N 



and likewise for S(P). This expression is consistent 
with the definition of a density operator S(r — ro)S(r,p) 
which, when transformed into the rotated frame yields 
the given result. Furthermore, all the given expressions 
can still explicitly depend on energy via the gauge in- 
variant parameter E, which however does not affect the 
discussion here. 



B. Operator equations of motion 



In order to calculate the densities (80) and (81), we 



need to evaluate the equations of motion for kinetic posi- 



tion R and momentum P. Using this result and (47), we 



I 

can now immediately write down the equations of motion 



dR 

~dt 



1 



[ST*] 



dP 



(1 + Gf) - s 



( P ) 



— = -m*p] = -(i 

dt in 1 1 1 u 



or, to leading order in fy. 
dR 



= 'dR 



dU 
dR 

dU 
f dP 



dt 

(1 - 0*") 



dp 



dt 



&H _ 
dP _ 
dH 
~dR 



dP 

lit 
dR 



(89) 

xS w , 
(90) 

(91) 
(92) 
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where we see that the effect of is related to a change 
of phase-space in the course of the diagonalization trans- 
formation. The beauty of this result is the symmetry 
in which effective magnetic and electric fields appear in 
these equations. 

We are now performing a quantum average of the ki- 
netic equations (91) and (92) with respect to a density 
matrix that is peaked around a certain value P c and R c , 
which in effect corresponds to a Gaussian wave-packet 
that is well beyond the limits of the Heisenberg uncer- 
tainty. Then our equations of motion (91) and (92) es- 



sentially become classical ones with all kinetic variables 
replaced by P c and R c . These equations of motion for the 
center of mass coordinates of a wave packet have been di- 
rectly obtained by Sundaram and coworkers^. The cir- 



cular current terms in (80) and (81) vanish for the wave 



packet, which is to be expected physically, since this de- 
scription reduces the electron to a point particle with 
coordinates P c and R c , and which does not posses any 
internal motion described by the circular currents. Math- 



ematically, after integration over the whole phase-space, 
one can rewrite them in terms of surface integrals which 
vanish for the well-localized wave-packet. 



C. Non-equilibrium description 

Now, it is just a matter of finding the density ma- 
trix p(R, P, t) in order to explicitly calculate anything 
within this framework, in particular, when one is inter- 
ested in non-equilibrium phenomena. One way to pro- 
ceed is within Keldysh formalism and to consider the 
kinetic equation for the lesser Greens function^, 



(93) 



which, after transformation into the rotated frame, can 
be evaluated to leading order correction in ^ by using 
( 47 1 which by virtue of ( 89 ) and ( 90 1 can be recast into 



di?, 

~dT 



dP, 

dt 



m 

~dt 



dP 



-IK 



o . 



(94) 



When we integrate the lesser Greens function over en- 
ergy, we obtain the density matrix 



p(R, P,t) 



dE 
2^ 



(95) 



and, since the terms proportional to 8eQ < vanish after 



integrating (94) over all energies, we are left with the 



following differential equation for the density matrix, 



dR 

~dT 



dP. 

~dt 



V P )p(R,P,t)=Q 



(96) 



which is nothing but a Boltzmann equation including 
quantum corrections in terms of Berry curvatures. In 
the end, it is quite analogous to the quantum Boltzmann 
equation along with the gradient expansion we used Ref. 
|4"2"1 Of course, the major difference is the absence of the 
collision integral in the present formulation, however it is 
possible to include it here as well, which is however left 
for future investigations. Of course, p can additionally in- 
clude discrete quantum degrees of freedom like the spin, 
which means that our rotated frame is block-diagonal and 
so is p. 

We now show how to obtain the distribution function, 
when it is known at some initial time, p(R, P, t = 0). In 
general, these problems can be formulated in terms of an 
inhomogeneous Boltzmann equation with a source term 



S(R,P,t), 



d t + — V fl 
dt 



dP. 

~dt 



Vp)p(R,P,t) = S(R,P,t) 



(97) 



To proceed in the usual way, we define the Greens func- 
tion g c , 

d t + — v R + — v P Jg c {R,p^R\p\t') 

= 5(t - t')S(R - R')5(P - P') , (98) 
for which we find the general solution 

Q c {R,P,t,R\P',t') 

= 0(t - 5{R - R c {t)) HP - Pc(t)) , (99) 

where O is the Heaviside step function, and R c (t) and 
P c (t) describe the classical orbits which obey the equa- 
tions of motion 



R c 

Pc 



dR 

lit 
dP 

~dT 



flc(0 = R' 

Pc(0 = P' 



This solution describes ballistic trajectories given by the 
kinetic equations ( 89 ) and ( 90 1 and which including the 



effects of the Berry curvatures, while the particle being 



15 



initially at phase space coordinates (R! , P' ,t'). For ex- 
ample, in the absence of any spatially dependent poten- 
tial, one simplyhas straight lines described by the Greens 
function Q c = Q(t-t ) S(R-V(t-t )) S(P-P'), where 



V = jjp- is the group velocity. This is analogous to 
other quasiclassical equations like the Eilenberger equa- 
tion which can be also described in terms of classical 
trajectories^ 21 . 

Now, the solution is readily given by 



/ + OO !• 
dt' J d d R' J d d P' g c (R,P,t,R!,P',t') S{R',P',t') , 



(100) 



where, due to the delta functions in Q c , we get con- 
tributions only from trajectories which end at (R, P) 
at time t. Since the trajectory is well defined by the 
given pair (R, P) already, we can find all corresponding 
points (R , P ) by going back in time so we can eventu- 
ally rewrite the above integral as 



p(R,P,t) 



dt' S(R c (t'),P c (t'),t') , (101) 



where now, the classical orbits are defined such that its 
position in phase space at time t is R c {t) — R and 

Pc(t)=P. 

Let us now go back to the original initial value prob- 
lem which we can easily solve by using the source term 
6(t)p(R, P,t — 0), so that the general solution reads 



p(R,p,t) = e(t) p(R' c (o),Km,o) 



(102) 



which has the desired properties for t > 0. This indeed 
corresponds to a motion in phase space as an incompress- 
ible fluid, as stated by Liouville's theorem. 



D. Polarization 

Another way to interpret the kinetic variables is ob- 
tained by considering the electrical polarization which 
we pursue according to the pioneering work of Vanderbilt 
and coworkers^. They derive the polarization in terms of 
adiabatic transport, where some perturbation that leads 
to a polarization in the crystal is adiabatically turned on. 
The current due to adiabatic transport is^ 



d d r d d p dR 
{2nh) d ™^lt 



BZ 



d d p 



Tr N £ 



VP) 



(103) 



where integration is only within the Brillouin zone. The 
electric polarization is then obtained by integrating over 
time, 



V(R) = J dtj(R,t) 



(104) 



The definition as an adiabatic transport process is nec- 
essary in order to obtain a truly gauge invariant results. 



In the periodic gauge howevei^, it is nevertheless possible 
to construct 



(2TTfl) d P ' 



(105) 



but in the general case, it is required to consider the full 
Berry connection structure (including in particular At) 
in order to obtain the correct result for the polarization. 
In the same spirit, one could define a polarization in mo- 
mentum space which can be thought of as a Dopplcr shift 
in the rotated frame or also termed anomalous velocity. 

In this sense, one could loosely interpret these Berry 
connections as shifts that the canonical variables acquire 
and which depend on full phase space, albeit these shifts 
are not directly physically observable, only when one does 
an integration with respect to either the position or the 
momentum variable one obtains the observable electric 
polarization in real space or reciprocal space, respec- 
tively. This is analogous to the quasi-probability dis- 
tribution p{r,p) in the Wigner picture, which can be 
given only physical interpretation as probability density 
when integrated either over whole momentum space or 
real space. However, if we try to transfer this idea, we 
still will have a shift that depends on the frame, i.e. it is 
still not gauge invariant. Only for a very specific gauge, 
one obtains the physical polarization. 



E. Bulk-boundary correspondence 



It is interesting to note that the circular current term 



d d P 



Tr N pW 



(106) 



appearing in the current density (87) is a manifestation 
of the bulk-boundary correspondence!^ In the homoge- 
neous bulk, these circular currents compensate each other 
and yield zero, but at the edge of the system, there will 
be residual currents flowing on the surface. Due to the 
circular nature of these edge currents, the associated edge 
states are also of topological nature. Therefore, a non- 
trivial (bulk) value of fi^ indicates a topological state of 
matter. 
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V. HIERARCHY OF EFFECTIVE THEORIES 
INCLUSION OF THE ELECTROMAGNETIC 
FIELD 



In many cases, one starts from an effective theory 
which requires dealing with gauge invariant momentum 
and/or position operators in the original frame. This 
might happen if one uses an effective theory derived 
from a more comprehensive theory; or more commonly, 
one simply wishes to include an external electromagnetic 
field. Qualitatively, one expects no new concepts or phe- 
nomena to emerge in effective theories at a lower level 
in the hierarchy, as one could also go directly from the 
topmost theory to the lowest in a single step. 

Here, our focus lies on the combined effect of the elec- 
tromagnetic field and the Berry connections emerging 
from the diagonalization of a band Hamiltonian. In this 
case, we are starting from a theory which is gauge invari- 
ant in real space, thus formulated in terms of the mini- 
mally coupled momentum p— qA{x). Since all our deriva- 
tions so far are only linear in h, we treat the electromag- 
netic vector potential qA (q is the charge of the particle, 
so we usually have q — — e) on a different level than HA, 
thus formally taking h and q as independent expansion 
parameters. To begin with, we have the usual canonical 
pair x and it and the kinetic pair x and Tr+qA(x), where 
A = (—</>, A r ) is the usual 4-component vector potential 
of the electromagnetic field. In the following, we denote 
B(r,t) — V r xA r and E(r,t) — -d t A r + V r A Q (not to 
be confused with the energy E) as the external electric 
and magnetic fields. 



The Hamiltonian is specified as H(x,p — qA r (x)) and 
is diagonalized by an appropriate unitary matrixW(a;,p— 
qA r (x)) 1 so that in the diagonalized frame we arrive at 
the kinetic pair of variables 

X = x — HA n (x,p — qA r ) 

U = TT + qA(X) + HA x (x,p-qA r ) . (107) 

Since derivatives now also act on A r , we have 

A Xk = A Xk - q(d Xk Ai)A Pl (108) 

so that we can write more explicitly 

R = r + hA P (109) 
E = e — q(f) + HAr + hqE ■ Ap (110) 
P = p - qA r (x) - KAr - hqB x A P . (Ill) 

Note that Ax, as well as Ap still carry the full gauge 
invariant momentum p—qA r and for simplicity, all Berry 
connections involved are assumed to be Abelian. 

The fictitious fields in the presence of an electromag- 
netic field can be derived by plugging these kinetic vari- 
ables into the commutator relations ( 32 1 and ( 43 1 and 



keeping all terms up to order 0(h) and quadratic in 
the fields 0(q 2 ). In order to obtain manifest gauge in- 
variant results, we express everything in terms of ki- 
netic variables by using d Xj = dxj — q{dx j Ai)dp l and 
E(X) = E(x) + (Apd r )E(x). Denoting primed quanti- 
ties as the Berry curvatures describing the combined ef- 
fect of external electromagnetic and the diagonalization, 
we can summarize the results as 



B ,(r) = qB(R) + B M + q [Ti^efO-B - B<cT} + q 2 (B^-B) B (112) 

B' (p) = S (p) (113) 

S' M =qE(X)+S M +q<cT p E + qBxS (p) + q 2 Bx(ExB w ) (114) 

£/(p) = £(p) +q ExB<- p) (115) 

g-7 = &% + q[(B™-B)8 ij -Bf>B j ] . (116) 



Here, Unprimed quantities B, S, correspond to the 
Berry curvatures of the system in the diagonalized frame 
in absence of the external electromagnetic field qA. 

It is important to note that these results have to be 
substituted into the equations of motion (89), ((90l and 



not (91 ), (92) since the latter has been obtained by drop- 



ping terms beyond leading order corrections in Berry cur- 
vatures, whereas our results for the fictitious fields are 
in fact higher order in Berry connections (q 2 h actually 
already corresponds to terms of third order in the expan- 
sion of the Moyal product). 

The term E x 2? (p) of £' {v) is the anomalous velocity 
term that is for example responsible for the quantum Hall 



effect, since the Berry curvature £? (p) becomes non-trivial 
when the system is in the Quantum Hall state. Similarly, 
in systems with spin-orbit interactions, this term consti- 
tutes the intrinsic contribution to the anomalous Hall 
effect. The reciprocal effect thereof is described by the 
term qJ3xf (p) , where a external magnetic field trans- 
forms a momentum space electric field into a real-space 
one. 

,7p gets modified by a magnetic field term that resem- 
bles a dipole interaction between a real-space and the mo- 
mentum space magnetic field. In fact, as we have seen 
in result (78| ), and as we will see explicitly at the end 
of this section, Q ,rp is directly related to the magnetic 
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dipole energy in the case of a two-level model, and thus, 
it is the magnetic field term in Q ,Tp which will give rise 
to the Zeeman energy in the Hamiltonian. Furthermore, 
the term quadratic in the fields can be absorbed by 9 ,7p , 



i.e. fictitious electric field can be recast into 



£' M = qE + £ M + qS'^E + qBx 



(117) 



From this expression, we see that the external electric 
field enters effectively as q(t + & TV )E, and the electric 
field acting on the system is thus renormalized. An anal- 
ogous renormalization of the external magnetic field ap- 
pears also in the result for . 



The modification of the transformation rule (481 for 



arbitrary observables is easily found by virtue of ( 108 1 



O' = O+[y[-y 1> 6 \+pE{A P ,d E O }-e ijk qB k ^{A Pi) [d Pj -^A Pj ,O ]} , 



(118) 



where again, primed symbols represent quantities in pres- 
ence of an external magnetic field. In particular, eval- 
uating the above expression for the Hamiltonian, we 
find for the magnetizat ion energy in the rotated frame 
H M = TrO' 7 " (sec Eq. (Tnl)), where 



n'.T = nr. 



q[{W-B)5 ij -£lYB j ] 



(119) 



Notice the resemblance between this expression and 



(116 1; in fact, it turns out we obtain £2' by substituting 



the Ber ry cu r vatur es with the analogue quantities of £2 in 
results ( 112 ]-( 116 ). In two dimensions, is zerop^and 
qB-3 ip) , i.e. a factor 2 less than 
Tre 7 ". Nev- 



thus Trg 7 ^ =~Tt£r p 

in the 3D case, and likewise, in ID TrQ ,rv 
ertheless, we see that the magnetic field couples in the 
shape of a dipole term _Bi£? (p) to the momentum space 
magnetic field 2? (p> . Since the magnetic dipole energy is 
already of order fc, we can equally write B(X) and B(x), 
since substituting canonical and kinetic variables affects 
only higher orders in fy. 



V 



where m is the rest mass and p the momentum of the 
electron. As 



a 



CT 

a 



= 



(121) 



H acts on the 4-dimensional Dirac Spinor and thus can 
be considered as a 4-band model. 

The Foldy-Wouthuysen transformational which brings 
the Dirac equation into diagonal form is described by the 
unitary matrix 



Wo 



(E P + mc 2 )l 4 + c/3aP 
^E P (E P + mc 2 ) 



(122) 



where Ep — c\J P 2 + m 2 c 2 is the relativistic energy of the 
electron with gauge invariant momentum P. Performing, 
as outlined before (cf. Eq. ( 22 1), the diagonalization with 
respect to Uq, we arrive at 



VI. THE DIRAC EQUATION - AN 
ALTERNATIVE PERSPECTIVE ON DERIVING 
THE PAULI-HAMILTONIAN 

The goal of this section is illustrate our formalism, 
while at the same time serving as a concise summary of 
the framework, essentially making use of all the central 
results. The Dirac equation and its descendants like the 
Pauli Hamiltonian are the most fundamental equations in 
condensed matter physics, nevertheless, we present some 
interesting insight into the physics revealed by this rather 
simple equation. 

We consider the Dirac Hamiltonian in the presence of 
a scalar potential and minimally coupled to the magnetic 
fielcP, 



H = ca(p - qA) + mc 2 /3 + V(x)l 4 



(120) 



H=U*H*U^ = E P a T z 



+ ( V(r)a + h{d r V)A { p ] + qEpB-B 



2<P) 



(123) 



The 2 x 2-matrices Tj denote the Pauli matrices in 
electron-positron space, while Ui describes the usual spin 
degree of fr eedo m. The magnetic dipole term is given ac- 
cording to (1191 in terms of an interaction term between 
magnetic field B = VxA and the fictitious momen- 
tum space magnetic field B lP \ to be specified below. We 
note that in this case, all the correction terms are equal 
for both electron and positron bands, albeit the term 
positron becomes only meaningful when all negative en- 
ergy states are completely occupied. 

The full matrix structure of the Berry connection, split 
into diagonal and off-diagonal parts, i.e. Ap = Ap tq + 
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A 



reads 



( d ) = C 2 PX(T 

p 2Ep(E P + rr 



X c Pxi 

Ah mc 



A 



(°) 



P(Per) c 3 



1Ep IE p(Ep 



2h 



(124) 
(125) 



where in the last step, we are dropping terms of order 
0(p/mc) 2 . 

The kinetic variables for both positive and negative 
energy states read explicitly 



P 
R 



P 
r 



qA(r) 



qhBxA ( p ] 



hA 



(d) 



so that the Berry curvature (34 1 is also equal for both 
bands, and yields 



„, , kmc 4 

B (p) = =-<r - 

2E% 



P(P-o-) he 4 

2Ep{E f 



mcr 



\2 
— — ^-<7 

2h 



(126) 



where again, in the last step we took the non-relativistic 
limit. Of course, the general expressions are still valid 
for arbitrary velocities, for example, in the opposite, the 
ultra-rclativistic limit, we find 



B 



( p) p>™ P(P-er)h 

2P 4 



Using the results (112 1 - (116), we immediately find 
f" = = and 



B (r) = qB(R) + q 2 (B™-B) B 

Q? ] =q[(B M -B)5 l] -B? ) B j \ 



(127) 



where here, we used E(R) = 0, since we are already 
using the scalar potential V(R). 

In terms of kinetic variables, the Hamiltonian T-L can 
now be rewritten (let us restrict ourselves to the positive 
energy branch, whose excitations correspond to electrons, 
so that q = — e), neglecting terms of order pi 2 , 



H = E P + V(R) + qE P B(R)-B 



(p) 



(128) 



so that in view of this, the spin-orbit interaction that ap- 
pears in the rotated frame can be reinterpreted as being 
a result of the shift the kinetic position operator attains. 
However, we should keep in mind that this shift is gauge 
dependent and becomes physically meaningful only when 
integrated over the whole momentum space (which then 
is equivalent to the polarization) or in the form of the 
Berry curvature, viz the momentum space magnetic field 
appearing in the kinetic equations. In the non-relativistic 
limit, the magnetic dipole energy is just the usual Zee- 
man term (using (71) and (72)) 



-H a B , 



(129) 



where we introduced the magnetic moment of the elec- 
tron spin, fi s = —g s ^cr/2 = — /ib"", and we assume 



a g-factor of 2 within the validity of the Dirac theory 
without quantum corrections from the radiative field, and 



/iB = ^ is the Bohr magneton. 



Using results (89) and (90) together with (127) and 



( 128 1, we obtain the operator equations of motion of a rel- 



ativistic electron which moves in an electromagnetic field 
that is smooth on the scale of the Compton wavelength 
A c . We do not explicitly write the expression here, how- 
ever, our result is consistent with the work of BliokrPH. 

For the sake of our discussion, we will nevertheless 
specify the equations of motion in absence of a magnetic 
field, 



dR _ dH 
~dt ~dP~ 

dp _ _m 

~dt ~ dR 



dt -dP + 2hdt XfT > 

= —eE(R) , (131) 



whose anomalous velocity term is a factor 2 larger than 
what one would obtain from the Hamiltonian (123), 



naively taking the canonical variable r to be the phys- 
ical position operator. Essentially, the same conclusion 
has been reached in Ref. [3] for the case of the relativistic 
Dirac equation. 



Th e particle density is, according to results (86) and 
p7], given by 



n(R,t) 



d 3 P 
{2irhY 



Tr 4 p(l - 2qB-B M ) , (132) 



and let us take the vacuum state with all negative en- 
ergies filled, so that p = Q(—r z Ep) and the first term 
gives the vacuum charge. However, from the correc- 
tion factor proportional to the magnetic field one might 
conclude that a sufficiently strong magnetic field leads 
to a charge accumulation with respect to the vacuum 
state. This kind of charge accumulation is also present 
in topological insulators^. Here however, factors like 
D = 1 - Tr6 rp = 1 - are relevant only for enor- 

mous magnetic fields of w 10 10 T, so we can safely drop 
all corrections due to G^. 

It is also insightful to study the current density in the 
presence of the electromagnetic field, because it gives a 
direct meaning to the divergence-free current terms in 
( 87 1 . For the simplicity of our discussion, we restrict our- 



K 



selves to the non-relativistic limit, where B (v) 

and Ep = mc 2 +P 2 /2m, and furthermore, we treat 
external fields E and B only to linear order. Since we 
want to study excitations of positive energy, we ignore 
the completely filled lower band. 
Then, our charge density is 



n(R,t) 



d d P 

(2nh) c 



so that the charge current density becomes with the help 
of ( 87 ) and ( 89 ) as well as the definition of the spin mag- 
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netic moment (c.f. Eq. (129)) 



j{R, t) 



d d P p 



(2irH) d m mc- 



: n(R,t) Ex/j, s 



H s xV R n(R,t) , (133) 



where the surface term vanishes since we assume our elec- 
tron momentum to be peaked around some mean value 
P c . The density p can be explicitly determined from 
equation (96 1, here however, we only integrate this equa- 
tion over all momenta and arrive at the continuity equa- 
tion 



d t n(R,t)+V R -j(R,t)=0 



(134) 



The first term in ( 133 ) is the usual definition of a non- 



relativistic current, and the second term is the anomalous 
velocity contribution due to the electric field. Like for the 
electron wave-packet considered in equation (130), using 



the non-relativistic Pauli Hamiltonian without the dis- 
tinction between canonical and kinetic variables would 
yield an anomalous velocity term that is just half of the 
correct value. However, the last term seems unusual at 
first as it constitutes a persistent current which itself 
gives rise to a magnetic moment, 

i j d'RRxj prir = -i J d'RRx(p,xV,n) 

where Af = 1 is the number of electrons. This result 
suggests that the notion of electron spin and this internal 
persistent current are just different interpretations of the 
same effect. For example, a magnetic field couples to this 
circular current via the magnetic moment it generates, 
and gives rise to the Zeeman energy (129). This comes 



close to the original suggestion of an internal rotation of 
the electron by Uhlenbeck and Goudsmit^, yet unlike 
their idea, it is not the motion of a solid object, instead 
it is a genuine quantum phenomena, and thus does not 
suffer the same deficiencies concerning rotation velocities 
of the electron that would have to be faster than the 
speed of light. 

As discussed before, the spin-orbit term in the Pauli 
Hamiltonian is not gauge invariant, it gives the correct 
energy spectrum, but not the correct equations of mo- 
tion. That is why it has not been realized for a long 
time, as the Pauli Hamiltonian has been mainly put to 
test with electronic spectra in atoms or solid matter. 
To explicitly show this point, we go into a different ro- 
tated frame so that the Berry connection is changed to 
A p — > A p + fiVpx(p), and we get an additional term in 
the Pauli Hamiltonian h(V p x)(d r V). Let us assume for 
the moment that we are working in momentum space and 
in the operator representation of quantum mechanics, 
where the Hamiltonian formally looks the same. Now, 
performing a local gauge transformation in momentum 



space by adding the phase factor e lx ^ to the wave func- 
tion, we obtain e'^V (r)^^ = V(r - hV P x) = 
V(r) — fi{Vpx){d r V) + 0{fc 2 ) which exactly cancels the 
additional term, and we are left with the original Hamil- 
tonian. Usually, one enforces gauge-invariance in real 
space to obtain the electromagnetic field. This is a an 
example of gauge invariance in reciprocal space which 
leads to the phenomena of spin-orbit interaction, or the 
anomalous velocity. 

In the above analysis, we needed to drop terms of order 
h 2 , though here, H is just a formal expansion parameter, 
the real relevant scale being the Compton wavelength 
A r = — . For example, here, we find that the order 

L mc * 1 3 

of HA^ is 4T ^f e2 = A c ^ which is to be compared with 
our position coordinate, where one realizes that correc- 
tions are indeed small unless the potential V(r) varies 
strongly on the scale of the Compton wavelength A c . 
Consequently, the Hamiltonian (123) still describes the 



correct quantum behavior for sufficiently smooth poten- 
tials and thus, the term semiclassical expansion seems 
misleading. However, for the case of the Coulomb delta- 
like potential of the nucleus of an atom it is obviously no 
longer true, and this leads to corrections like the Darwin 
term which is essential to understand atomic spectrEp^. 

Sometimes, one interprets the terms appearing in the 
Pauli Hamiltonian as being a result of what one termed 
Zitterbewegung, and which one envisages as helical mo- 
tion of the electron. Then the spin-orbit interaction and 
its physical consequences like the anomalous velocity can 
be understood as a result of this motion. Obviously, we 
do not observe such an oscillatory motion here in the 
effective theory, when the excitation is confined to a sin- 
gle band. Essentially, due to the Wigner transformation, 
this rapid oscillation has been transformed into energy 
space, and instead, the effects like anomalous velocity 
now appear in form of inter-band corrections. We re- 
mark that in the case of massless Dirac-Fermions and 
if one considers excitations in the vicinity of the Dirac 
point, both positive and negative energy states have to 
be taken into account simultaneously so that excitations 
are not confined to positive or negative energies alone. 
In this scenario as for example in Graphene, one obtains 
those features termed Zitterbewegunap^. 



VII. CONCLUSIONS AND OUTLOOK 

In this work, we first studied a very general Hamilto- 
nian that contains an additional matrix structure describ- 
ing different bands, and the goal was to bring this Hamil- 
tonian into band-diagonal form which has been achieved 
by performing a rotation in band-space. The diagonal 
representation is very practical when one wants to study 
the low-energy response of the system because then, only 
one or few degenerate bands are relevant. Since in the 
very general case this diagonalization involves the pair 
of non-commutating position r and momentum p oper- 
ators, we performed the diagonalization perturbatively 
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in h by using the Wigner representation. We investi- 
gated the Hamiltonian and physical observables in the 
rotated frame and how Berry curvatures emerge natu- 
rally during this diagonalization procedure. Essentially, 
Berry connections describe corrections to the canonical 
variables like position and momentum due to inter-band 
scattering. This led to the distinction between canonical 
and kinetic variables, where the kinetic variables in the 
rotated frame are strongly linked to the canonical vari- 
ables of the original theory. The canonical ones are, on 
the other side, responsible for the proper quantum struc- 
ture due to their canonical commutation relations which 
defines the quantization. Therefore, both the canonical 
and kinetic variables are an important part of our effec- 
tive description, and only then will it be consistent. We 
also established a link between gauge invariance in mo- 
mentum space and the spin-orbit interaction, so that the 
Hamiltonian expressed in terms of canonical variables is 
enough to study energy spectrum of the system. 

Having in mind a gauge-invariant description, we ex- 
pressed the Hamiltonian and observables in terms of ki- 
netic variables, which naturally leads to the appearance 
of Berry curvatures which are gauge invariant and de- 
scribe various effects intrinsic to the band structure. In 
addition to Berry curvatures, we identified further gauge 
invariant objects which are related to circular, or persis- 
tent currents, which itself interact with magnetic fields 
or with each other and give for example rise to energy 
terms, like the Zeeman interaction. For the Dirac equa- 
tion, the electron is naturally delocalized in space and 
we found the appearance of an internal motion in the 
form of circular persistent currents giving rise to a mag- 
netic dipole moment, and which exactly corresponds to 
the spin magnetic moment of the electron. 

Using these kinetic variables, we formulated a Boltz- 
mann transport equation that incorporates intrinsic ef- 
fects expressed in terms of Berry curvatures. To this end, 
we consistently defined various quasi-probability densi- 
ties which are connected by a conservation law: Li- 
ouville's theorem in phase space. These can be used 
to obtain physically meaningful densities and currents 
within the effective description. Furthermore, it is rather 
straightforward to include impurity scattering in the 
same manner as in Ref. |42j which is subject of future 
work. The results derived in this work provide a good 
starting point to transform the collision integral into the 
diagonalized frame, automatically incorporating the ef- 
fect of an external electromagnetic field. 

A related work by Wong and Tserkovnyak treated a 
general 2-band system by studying a quantum kinetic 
equation approach in the rotated framed However, the 
rotation diagonalizes the Hamiltonian only to zeroth or- 
der in a gradient expansion, so that when studying gra- 
dient corrections the diagonalization is no longer exact. 
It is therefore no longer sufficient to consider only the 
diagonal part of the kinetic equation. To this end, it is 
necessary to take into account the off-diagonal elements 
which is in contrast to our formulation where already the 



rotation transformation includes gradient corrections. 

Despite performing a formal expansion in h, similar 
to treating the semiclassical limit, we never abandon the 
quantum description. One major advantage compared to 
other semiclassical approaches is that there is no need for 
a re-quantization (see Appendix [B] for some details). In 
fact, our real expansion parameter might be a different 
one, like the Compton wave-length as we have seen in the 
case of the relativistic Dirac equation. Also, our approach 
is systematic in the sense that we can go to arbitrary 
order in h or inter-band coupling. In this respect, it is also 
very interesting to look at terms second order in h which 
yield important contributions, for example the Darwin 
term in the case of the low-energy limit of the relativistic 
Dirac equation. Furthermore, new physical phenomena 
emerge at 0(h 2 ), like the magneto-electric coupling in 
insulators which received new attention r ecent ly, also due 
to the discovery of topological insulator j 51 * 52 l 

Then, one could study a variety of systems that involve 
both spin-orbit interaction (SOI) and an inhomogeneous 
and time-dependent magnetization. This is interesting 
for spin-orbit coupled semiconductors where one has dif- 
ferent types like Rashba SOI, Dresselhaus SOI or in the 
case of strong SIO in III-V ferromagnetic semiconductors 
one can utilize the Luttinger Hamiltonian for hole trans- 
port. In these systems the electron or hole spin is in- 
tricately coupled to both momentum and orbital degrees 
of freedom. In addition, one can have non-magnetic and 
magnetic impurity scattering that adds additional com- 
plexity to the system and is naturally studied in terms 
of a collision integral within the Boltzmann approach. 
Thus, the formalism developed in this work seems prac- 
tical to attach these kinds of systems in order to study 
various transport and dynamical properties. Finally, it 
would be interesting to study many particle effects like 
electron-electron interations. 

This work was financially supported by the DFG 
though SFB 767 and SP1285. 



Appendix A: Wigner Representation 

In the operator representation of quantum mechanics, 
we can write the canonical commutation relation (quanti- 
ties with hats indicate operators in the original picture), 



[Ttn,x v ] = ihS " . 



(Al) 



In 4-component notation, the Wigner transformation 
explicitly reads 

J d 4 z e™l % A (x + ~, x - |) = A(x, tt) , (A2) 

and the Moyal product obeys the axiom of associativity, 
A* B * C = A* (B * C) = (A* B) * C , (A3) 
and also 

(A * B) t = B f * At , (A4) 
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however, the Moyal product is non-commutative, i.e. A* 
B B * A. In fact, the latter property encodes the non- 
commutativity of operators in quantum theory which one 
can explicitly see by considering the commutator of two 
general observables dependent on x and 7r and which are 
band-diagonal, 

[A(x,ir) * B(x,ir)] = 2iA(x,n) sin(M/2) B(x,ir) 

= ih{A,B} p + 0(h 3 ) . (A5) 

It becomes essentially an extended Poisson bracket in the 
semiclassical limit which reduces to the normal Poisson 
bracket when neither A nor B explicitly depend on e 
which is usually the case. 

The Wigner representation has the very practical fea- 
ture of automatically symmetrizing operators, when one 
transforms back into the operator formulation, for exam- 



ple using the back transformation of ( A2 ) 



A{x- 



one obtains 



d 4 7T 



-iZTT jh 



A(x,tt) , (A6) 



xp —> — {xp + px) 



(A7) 



so what is usually required to be put in by hand in usual 
quantum mechanics in order to obtain Hermitian observ- 
ables, comes out automatically. Another, more direct 
way to see this is by using identities of the form 



(x * p + p * x) = xp . 



(A8) 



The special operator ordering obtained in this manner is 
the so-called Wigner- Weyl ordering. 



Appendix B: Some Notes on Requantization 

In this section, we address the question of the quan- 
tum theoretical aspects of our semiclassical analysis. In 
particular, in literature there is the question of how to 
properly quantize the theory in the Lagrangian formal- 
ism, since one needs to identify the canonical pair of vari- 
ables which is the starting point of the quantization. As 
mentioned previously, using wave packet analysis^, one 
obtains the same set of equations of motion for the cen- 
ter of mass coordinates of the wave-packet as in (91 ) and 



( 92 1 . Then one can find the Lagrangian, but one does 



not know anything about canonical variables, and in the 
general case, it is not always easy or possible to find the 
pair of canonical variables^. 

However, in our framework, we did not encounter such 
problems, since firstly, we are working in the Hamilton 
formalism, where one knows the canonical variables and 
secondly, our theory is quantized at any time, albeit ne- 
glecting terms of order 0(fi?). This quantization is en- 
coded in the structure of the Moyal product, or more 
explicitly in the special structure of the commutator in 
(47). Furthermore, without major effort, we can go back 



to the operator representation of quantum mechanics by 
undoing the Wigner transformation. 



Since many of our expressions (for example (19), (48), 
(29), etc.) are straightforwardly transformed back into 



the operator formulation of quantum mechanics, we for- 
mally obtain the same expression, except mixtures of r 
and p appear properly symmetrized (see Appendix [A]). In 
fact, we retain the full quantum theory which gives cor- 
rect results at least to order H. One just has pay attention 
that the Wigner transformation is only with respect to 
the canonical pair of variables so kinetic variables have 
to be replaced accordingly. 



* christian .wickles@uni-konstanz.de 

t wolfgang. belzig@uni-konstanz. del 

1 M. V. Berry, Proceedings of the Roy al Society of London!] 

A. Mathematical and Physical Sciences 392, 45 (1984)| 
' G. P. M ikitik and Y. V. Sharlai, |Phys. Rev. Lett. 82, 2147| 



(1999). 



CP 



D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82 



1959 (2010) 

J. N. Fuc hs, F. Piechon, M. O. Goerbig, and G. Mon- 
tambaux, The European Physical Journal B - Condensed 



Matter and Complex Systems 77, 351 (2010) 



K. M. D. Hals, A. K. Nguyen, X. Waintal, and A. Brataas, 



Phys. Rev. Lett. 105, 207204 (2010) 



R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47. 



1651 (1993)| 

F. D. M. Haldane, |Phys. Rev. Lett. 93, 206602 (2004)1 
N. A. Sinitsyn, Journal of Physics: Condensed Matter 20, 



023201 (2008) 



a N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 

N. P. Ong, |Rev. Mod. Phys. 82 1539 (2010H 
10 M. N. L euenberger and D. Loss, |Phys. Rev. B 63, 054414| 
(2001) 

12 



G. Sundaram and Q. Niu, Phys. Rev. B 59, 14915 (1999) 



D. Culcer, Y. Yao, and Q. Niu, Phys. Rev. B 72. 085110 

I ( 2005 )l . 

Li P. Gosselin, A. Brard, and H. Mohrbach, The European 

Physical Journal B - Condensed Matter and Complex Sys- 



tems 58 
L * P. Gosselin 



137 (2007)| 

Hanssen, 



and H. Mohrbach. 



Phys. Rev. D 



77, 085008 (2008) 

P. Gosselin, H. Boumrar, and H . Mohrbach, EPL (Euro 
physics Letters) 84, 50002 (2008)1 
10 C. H. W ong and Y. Tserkovnyak, Ph ys. Rev. B 84, 115209 
(2011)1 



22 



17 D. Culcer and Q. Niu, |Phys. Rev. B 74, 035209 (2006) 

18 C. Gorini, P. Schwab, R. Raimondi, and A. L. Shelankov, 



19 



Phys. Rev. B 82, 195316 (2010) 



R. Winkler, Spin-orbit Coupling Effects in Two- 
dimensional Electron and Hole Systems (Springer- Verlag 
Berlin Heidelberg, 2003). 



20 E. Wigner, |Phys. Rev. 40, 749 (1932). 

21 G. Tat ara and H. Kohno, |Phys. Rev. Lett. 92, 086601 



(2004) 



G. Tatara, H. Kohno, and J. Shibata, Physics Reports 
468, 213 (2008)| 
^ C. H. W ong and Y. Tserkovnyak, |Phys. Rev. B 80, 184411 
(2009)| 



G. E. W. Bau er, E. Saitoh, and B. J. van Wees, Nature 
Materials 11, 391 (2012)] 
^ D. Hsieh, D. Qian L. Wray, Y. Xia, Y. S. H or, R. J. Cava, 
and M. Z. Hasan, Nature 452, 970 (2008) 



X.-L. Q i and S.-C. Zhang, |Rev. ModTP hys. 83, 1057 
(2011)| 



27 X.-L. Q i, T. L. Hughes, and S.-C. Zhang, Nat Phys 4, 273 
(2008) 

38 rp T7~T 



T. Yokoyama, J. Zang, and N. Nagaosa, Phys. Rev. B 81, 



241410 (2010) 

:v K. Nomura and N. Nagaosa, |Phys. Rev. B 82, 161401 



(2010) 



C. Wickles and W. Belzig, Phys. Rev. B 86, 035151 (2012) 



Y. Tserkovnyak and D. Loss, Phys. Rev. Lett. 108, 187201 



(2012)| 

^ M. Matsuo, J. Ieda, E. Saitoh, and S. Maekawa, jPhys 



Rev. Lett. 106, 076601 (2011) 
00 M. Matsuo, J. Ieda, E. S aitoh, and S. Maekawa, |Phys 
Rev. B 84, 104410 (2011) 

■» "V TV,™.!,™™,,™!, „«J /~< TJ 



Y. Tserkovnyak and C. H. Wong, Phys. Rev. B 79, 014402 
(2009) 



3b L. L. F oldy and S. A. Wouthuysen, |Phys. Rev. 78, 29 
(1950) 



JU J. Rammer and H. Smith, |Rev. Mod. Phys. 58, 323 (1986) 



37 J. E. Moyal, Mathematical Proceedings of the Cambridge 


Philosophical Society 45, 99 (1949) 




J ° S. E. Barnes and S. Maekawa, Phys. 


Rev. Lett. 98, 246601 



(2007)| 

3 " R. Shindou and L. Balents, |Phys. RevT 



B 77, 035110 



J. Shi, and Q. Niu, Phys. Rev. Lett. 95, 137204 



(2008) 
D. Xiao, 
(2005) 

Since effectively projects the Berry connection ma- 

trices _4 (d) , so that we have [S (d) ,.4 (d) ] = 0, the Taylor 
expansion is mathematically well-defined. In the case of 
Abelian gauge fields, the situation is trivial, while in the 
non-Abelian case, we need to make the assumption not to 
probe the internal, degenerate structure and accordingly, 
So should be diagonal in the corresponding sub-space in 
order to describe an observable of this kind. 



C. Wickles and W. Belzig, Phys. Rev. B 80, 104435 (2009) 



A. Shelankov, Journal of Low Temperature Physics 60, 29 
(1985)| 

** D. J. Thouless, |Phys. Rev. B 27, 6083 (1983)1 
45 A. M Essin and V. Gurarie, |Phys. Rev. B 84, 125132 

However, even in lower dimensional system we are still 
using all t hree d imensions in our vector calculus and ac- 
cording to (116 1, 7^ 0, so we have to pay attention to 
these subtleties for dimensions < 2. 

F. Schwabl, Advanced Quantum Mechanics, edited by 
R. Hilton and A. Lahee (Springer- Verlag Berlin Heidel- 
berg, 2008). 

K. Y. Bhokh, |EPL (Europhysics Letters) 72, 7 (2005)] 

G. E. Uhlenbec kand S. Goudsmit, Nature 117, 2 64 (1926). 
M. Katsnelson, |Eur. Phys. J. B 51, 157 (2006)| 
A. M. Essin, A. M. Turner, J. E. Moore, and D. Vander- 



M. Essin, A. M. Turner, J. E. Mo ore, 
bilt, Phys. Rev. B 81, 205104 (2010)| 
52 S. Coh, D. Vanderbilt, A. M alashevich, and I. Souza, |PhysT 
Rev. B 83, 085108 (2011)] 



